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The interaction between a superconducting vortex or antivortex in a superconducting film and a 
, magnetic dipole with in- or out-of-plane magnetization is investigated within the London approxima- 

tion. The dependence of the interaction energy on the dipole-vortex distance and the film thickness 
is studied and analytical results are obtained in limiting cases. We show how the short range interac- 
tion with the magnetic dipole makes the co-existence of vortices and antivortices possible. Different 
configurations with vortices and antivortices are investigated. 
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■ I. INTRODUCTION 

^3 ■ 

c/2 ' The value of the critical current is one of the decisive factors for the usefulness of a superconductor (SC) . Large values 
of the critical currents in superconductors are usually obtained through pinning of vortices to different inhomogeneities 
^ in a superconductor. In this respect, external and internal surfaces can be treated as inhomogeneities. □ Pinning centra 
^ , are introduced e.g. by locally destroying the crystal structure through e.g. bombardment with high energy particles. 
' , ' Recently, nanostructuring of a superconducting film with a regular area of holeso has lead to a large increase of the 
^ \ critical current, in particular at the so-called matching magnetic fields. An alternative approach is to deposit an 
Q t area of ferromagnetic dot&Jiear a superconducting film which acts as very effective trapping centra for the vortices. 
' Recently, it was predictedp although not verified, that an increase of the pinning effects by two orders of magiutude 
can be realized. After substantial progress in the preparation of regular magnetic arrays on superconductors □ and 
I ' considering the importance of such structures for magnetic device and storage technologies, these hybrid systems 
^ , became very interesting both from a theoretijcal and an experimental point of view. Macroscopic pinning phenomena 
' have already been explored experimentally, uU but a theoretical analysis of the piagnetic and superconducting response 
T-H \ in such systems is still lacking. In the majority of recently proposed modelallilil the inhomogeneous magnetic field of 

■ a ferromagnet induces screening currents in a SC, which, in turn, generate a magnetic field infiuencing the applied 
' field. Consequently, this problem must be solved self-consistently. Furthermore, the finite thickness of both the 

superconductor and the ferromagnetic material was not taken into account in previous theoretical treatments. 

Other theoretical studies involving finite size ferromagnets were maiuUj_c.estricted to the problem of a magnetic dot 
-j~>, : with out-of-plane magnetization embedded in a superconducting film.E3Ej Marmorkos et alL3 solved the non-linear 
' Ginzburg-Landau (GL) equation numerically, with appropriate boundary conditions for an infinitely long ferromag- 
\ netic cylinder penetrating the superconducting film, and found a correspondence between the value of the magnetiza- 
I 1 tion and the vorticity of the most energetically favorable giant-vortex state. At that time, the interaction between a 
' O ' small ferromagaejtic particle, which may be considered as a point magnetic dipole (MD), and a type II superconductor 
C : was of interest,oll3 and the magnetostatic energy ^md levitation force acting on a dipole were calculated, but still 
^ within a thin film approximation. In experiments,t3 the magnetic field of the MD was supposed to be weak, and 
. . not able to drastically change the structure of the superconducting state in the sample. Obviously, the MDj^ould 
create additional vortices near the surface, and this process can be described theoretically as proposed in Ref.cJ, by 
a simple comparison of the free energies of the system with and without vortex. However, the spontaneous creation 
of a vortex-antivortex pair as a possibleJ,ower energy state was never considered. 
^ Motivated by recent experiments, EJila we focus in this paper on a model system consisting of a type II thin 
superconducting film (SC) and a magnetic dipole placed above (below) it which acts as a pinning center (Fig. ^). We 
study in detail, using the London approximation, how the system is perturbed in the neighborhood of the dipole. The 
superconducting film lies in the z — Q xy-plane while the MD is positioned at {x,y,z) = (0,0, a), and is magnetized 
in the positive z- or x-axis direction. To avoid the proximity effect and exchange of electrons between MD and SC we 
suppose a thin layer of insulating oxide between them as is usually the case in the experiment. 

The paper is organized as follows. In the next section we present the general formalism. In Sec. Ill, we discuss the 
pinning potential of the magnetic dipole. In Sec. IV the total interaction energy in the system is calculated in the 
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presence of vortex-antivortex pairs, and their most favorable position is determined. The question of stabihty of such 
vortex configurations is then analyzed and our conclusions are given in Sec. V. 



II. THEORETICAL FORMALISM 



We consider a magnetic point dipole with magnetization M, placed outside a type II SC film interacting with a 
single vortex in the_SC. Within the London approximation, the total energy of the stationary magnet-superconductor 
system is given byt£l 
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where A is the penetration depth and h is the total field present in the system. The integral of the first term is taken 
over the volume of the SC, and the integral of the second term is taken outside, except for the volume of the dipole. 
Therefore, we may rewrite Eq. (|l|) as 
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where we integrate over the volume inside y*^'^ and outside V^°'^ the superconductor, while denotes the volume 

of the magnetic dipole. We choose the surface S far away from the superconductor where we can apply the boundary 
condition h 0. Due to the London equation, the field of the magnetic dipole satisfies 



hm + }?rot [ roth„i ) = 



(2) 



The magnetic field and corresponding vector potential can be written as h = hm + hy, A = Am + Ay where indexes v 
and m refer to the vortex and ferromagnet, respectively. Now, from Eq. we obtain the vortex-dipole interaction 
energy as 
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which after simple integral transformations becomes 
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where = ($p, = (0, $o/(27r/9), 0) denotes the vortex magnetic flux vector ($o is the flux quantum), and M 

is the magnetization of the dipole. As one can see, the interaction energy in this system consists of two parts, namely, 
(i) the interaction between the Meissner currents generated in the SC (jm) by the MD and the vortex, and (ii) the 
interaction between the vortex magnetic field and the MD. 

In order to obtain the value of the current induced in the superconductor by a mam£tic dipole with moment m 
(M — rh5{x)5{y)5{z — a)), we solve the corresponding equation for the vector potentiaO 



rot (^rotAm^ + -^8 {d/2 - \z\) Am = 'inrot {rhS{x)S{y)S{z - a)) , 



(4) 



where S{...) is the Dirac delta function. The results for both vector potential and magnetic field, for different orien- 
tations of the magnetic moment of the dipole, are given in the Appendix. 



III. MAGNETIC DIPOLE - VORTEX INTERACTION ENERGY 



The superconducting current induced in an infinite superconducting film with thickness d {—^ < z < ^) hy a 
magnetic dipole with out-of-plane magnetization^ i.e. rn = mcz located at z = a is obtained as a solution of Eq. (Q) 
which is given by 
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z) = -^^^ J dqexpi^-q(^\a\-^]'fq^Ji{qp)C{q,z), (5a) 




with 



C(^q z)^ coMHi + z))+q smh{k{^ + z)) ^^^^ 
' (fc^ + q^) sinh(kd) + 2fcg cosh{kd) ' 

where fc = -^/l + 9^, p = \/ +2/^, smft, and cos/i denote the hyperbohc trigonometric functions and Jv(x) is the 
Bessel function. For a MD placed under the SC, one should replace z by — z in Eq. (5b). The magnetic moment of 
the dipole is measured in units of mg = $0-^1 and all distances are scaled in units of A. These units will be used in 
the rest of the paper. Clearly, Eq. can also be used for a film of finite dimensions L in the limits p,a « L. The 
above integral (pa) can be solved analytically for certain asymptotic limits 
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, p >> max{l, l/d) or / >> max(p, l/d), 
p,l< i, d « 1, 



where I — \a\ ~ ^ is the distance between the MD and the top surface of the SC film. 

Inserting the well known expression for the magnetic field of a vortex (see Eq. (A4) in the Appendix) , into Eq. (|^) 
we find 

Fmv = ^Ui_ (p„) , (7a) 

where p — pv denotes the position of the vortex, and 
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In some limiting cases it is possible to solve this integral analytically: 

1 /2 

(1) for d << 1 and py < l/d we found: U± (py) « -d/{2 {pi + P) ' ); 

(2) for d < 1 and a) > a: Ui_ (p,) « -d/(2 {pi + I'^f'^) + i^d^ [Ho{pyd/2) - Yo{pyd/2)] /8, 

b) Py » l/d: Ui_ (py) « -2/{dpl)- 

(3) for d > 1 and py » max(^, l/d) or I » max(p^, l/d) we obtained: U±_ {py) Ki ~l j {pi + P^^^"^ . 

Here, Hyix) and Yy[x) denote the Struve and Bessel function, resjjwsptively. One should notice that our asymptotic 
results for a thin SC filtu differ from those given in Eq. (7) of Ref.til by a factor of two, and are in agreement with 
the calculations of Ref.y (Eq. (3.16) of Ref.Ea). 

In Fig. |^(a,b) the interaction energy as function of the position of the vortex are shown, for different (a) vertical 
positions of the magnetic dipole, and (b) thickness of the SC. Please note that the asymptotic expressions for an 
extremely thin SC, namely d « 1, give a very good description of the interaction potential (see inset of Fig. ||(b)). 
As expected, the problem of the interaction energy in this system is axially symmetric and this is illustrated in 
Fig. |^(c). In Fig. ^(a), we notice that the energy curves for different vertical positions of the dipole cross each other 
for large MD-vortex distances, for the case of a thick SC, namely d> 1. Consequently, for a MD which is higher above 
the SC film, interaction with the vortex is stronger at large distances, i.e. py » l/d (see our approximate results 
above, case (3)). This can be understood as follows. In the case of a thick superconductor, the magnetic field does 
not penetrate through the SC, and the component of the field tangential to the surface becomes important. When 
the MD approaches the SC surface, the value of this field component at large distances from the dipole increases, as 
shown in the inset of Fig. |^(a). Therefore, the interaction between the dipole and the vortex is stronger for smaller I 
when the MD-vortex distance py is small, and opposite, for large py, the interaction energy grows if I is larger. 

Naturally, the minimum of the interaction energy depicts the energetically favorable position of the vortex. There- 
fore, from Fig. |^ one notices that for parallel alignment of the MD magnetization and the vortex the interaction will 
be purely attractive, and the vortex is pinned under the MD, independently of the parameters of the system (thickness 
of the SC, vertical position of the MD and its magnetization strength, etc.). These parameters only determine the 
strength of this attractive interaction. For the anti-parallel orientation of the magnetic moment, the vortex is repelled. 
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Following the same procedure, for a MD with in-plane magnetization (m = mCp, for (p = 0, and m = for </? 7^ 0) 
we have 

where (p is the angle in the SC plane between the direction of the magnetic moment of the dipole and the point of 
interest, and 

00 

^11 ^J'^ h{k + rcoMkd/2)/ ' ('ZP.)exp(-<zO, (8b) 



for which we obtain the following analytic asymptotic results 
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, p » max(Z, l/d) or I » max(/9, l/d), d> 1, 
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In Fig. ^(a,b), the energy is plotted as function of the position of the vortex when it is moved along the direction 
of the MD magnetization ((^ = tt), and Fig. ||(c) shows the complete contour plot. The vortex is attracted to one 
side (where the magnetic field of the MD penetrates the SC in the same direction as the vortex, see Fig. 0(b)) and 
repelled on the opposite side of the MD. Moreover, it is pinned in a spot whose position does not depend on the 
magnetization of the MD but it is influenced by the vertical position of the dipole and the thickness of the SC. When 
the dipole approaches the SC or the thickness of the superconductor is increased, the interaction becomes stronger 
and the vortex will move closer to the MD. 

From our asymptotic expressions for the interaction energy the most favorable position of the vortex can be easily 
obtained by minimizing the energy. In the thin film approximation (d << 1), we obtain pi = ^2 + 2^5^/2 « 1.272 Z, 
while for d > 1, but I » l/d we have p* = l/\/2 « 0.7071 I. In Fig. |^(a) we show the d— dependence of the 
equilibrium position as function of the thickness of the SC. The p*(0 function is in general non-linear, as shown in 
Fig. ^(b). A peculiar fact is that for d > 2, the energetically preferable position of the vortex depends only on /. This 
follows from the fact that the magnetic field of the dipole penetrates only in the surface layer of the superconductor 
down to the penetration depth A so that further increase of the SC thickness does not affect the total interaction. 
From Fig. ^(b) one should notice that our asymptotic results describe nicely the behavior of the system for d << 1 
or I » d, when d > 1. 

For the case of in-plane MD, it should be emphasized that the vortex is attracted to the side of the MD where the 
magnetic field is parallel to the field of the vortex, which is similar to the case of the out-of-plane magnetized dipole. 
This conclusion also follows from the observation of the interaction between the induced currents and the vortex. For 
in-plane magnetized MD above the superconductor, the components of the superconducting current in the absence of 
vortices are given by 

jp{p, V^z)^- ^"^2^°^^^ ^ J dqexp{-ql)qJi{qp)C{q,z), (10a) 
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j,(p,(p,z) = 0. (10c) 



We obtained the following asymptotic behavior of these components: 
a) p » max(^, l/d) or I » max(p, l/d) 
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b) p,/ < i, d « 1 
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For the case of the MD under the SC, these currents change sign. The detailed expressions for the magnetic fields 
and the vector potential as solution of Eq. (^) are given in the Appendix. The vector plots of the screening currents 
for both directions of the MD jMagnetization are presented in Fig. ^j. 

Previous theoretical studieaHJ on the interaction of a magnetic moment and a vortex used a Gibbs free energy 
argument to show that this force was attractive. In our case, the pinning force between a vortex and the MD consists 
of two parts: (i) the interaction between the vortex fringing field and the magnetic moment, and (ii) the interaction 
between the vortex and the screening currents created in response to the magnetic dipole field. One should notice 
(Fig. ||(a)) that although the magnetic field of a MD with out-of plane magnetization changes sign in the SC plane, 
the superconducting screening current is always in the same direction, and therefore, the interaction between the MD 
and a vortex with parallel orientation of the field is attractive (consequently, the vortex will sit under the MD), and 
analogously, for anti-parallel orientation of the magnetic fields they repel each other. 

For the case of a MD with in-plane magnetization, the Meissner current vector plot in Fig. |5|(b) shows that the 
direction of the current is such that the vortex is attracted to the region where the magnetic field of a MD is parallel 
to it and repelled on the other side of a MD due to the opposite direction of the current. 



IV. PINNING OF MULTIPLE VORTICES AND/OR ANTI- VORTICES 

A. In-plane magnetized dipole 

Obviously, in the case of in-plane magnetization, the vortex will be repelled on the side of the MD to which the 
magnetic moment points to and pinned on the opposite side (and vice versa fop-the antivortex). The case when the 
dipole itself creates vortices was experimentally investigated by Van Bael et alS3, where it was found that the stray 
field of the in-plane magnetic dipole induces a vortex-antivortex pair (VA) at the poles, at positions predicted by 
our analysis (vortex towards negative pole, antivortex towards positive pole). Therefore, this asymmetric pinning 
potential provides stability for vortex-antivortex configurations. One would expect that the position of the pinning 
sites depends on the magnetic moment of the MD and the number of fiux quanta carried by each vortex. To determine 
the preferable position of the vortex-antivortex pair in the presence of an in-plane MD we put the vortex at {py, tpy) 
and the antivortex at (pau, Vav) with respect to the MD and calculate their interaction with the dipole. The vortex- 
antivortex interaction energy is given byo 

Fy^UR,Ly,L,y)^-^^^^^^ (^dKo{R) + 2 Jdql^^^ , (12) 

where R = \/ p1+ Pav ~ '2pyPavCos{ipy — ipav) is the distance between the vortex and antivortex, Kq (x) is the Mac- 
Donald function, and Ly, Lav are the vorticity of the vortex and anti- vortex, respectively. 

First, let us suppose that the vortex-antivortex pair appears due to the stray field of the dipole. In that case, 
this pair nucleates where the superconducting current is maximal, namely, under the dipole (Fig. |^(b)), with vortex 
and antivortex situated on opposite sides of the dipole. Therefore, due to the symmetry of the applied potential 
(Fig. |3|(c)), we may assume that py = pav, and ify = tt — ifav To investigate the stability of such a pair of vortices, 
we calculate the total interaction energy in this system with respect to the position of the vortices. The results are 
shown as contourplots in Figs, ^a-c) for different values of the magnetic moment of the dipole. The lowest value of 
the interaction energy gives the energetically favorable position of the vortex-antivortex pair [apart from the global 
minimum at a; = (zero distance between vortices), where the energy equals — oo]. Obviously, the lateral motion of 
vortices is forbidden by the high energy barriers, while along the direction of the dipole, also an energy barrier exists, 
but lower than the lateral ones, showing us the preferable direction of possible vortex-antivortex annihilation. 

Knowing this, we put the vortex and antivortex along the direction of the MD, in such a way that py = pay, and 
(fy = TT, ifav = 0, and investigate the pinning potential as function of the distance of the vortex and antivortex from 
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the MD. We substitute R = 2py in E n. ( [T^) and add this energy to the values of the MD- vortex and MD-antivortex 
attraction energy obtained from Eq. (|8a|). The result is shown in Fig. ^ for (a) different values of the vorticity {Ly, 
Lav), and (b) different magnetization and Ly = Lav = 1- After minimization of the total energy over py, we obtain the 
position of the vortex-antivortex pair. In the thin superconductor case (rf < 1) this leads to the non-linear equation 



2'Kmpy 




Ly + Lay 



(13) 



Although all forces acting in the system are attractive, it is clearly visible that annihilation of vortices can be 
prevented by their strong confinement at the poles of the MD (the position of vortex-antivortex pairs is illustrated by 
the open circles in Fig. [^). The interaction energy shows an absolute minimum for py = 0, implying the annihilation, 
but in order to do that vortices have to cross a potential barrier. Therefore, the vortex-antivortex configuration can 
be stable in this system. To determine if this vortex state is metastable or the ground state one should compare its 
total energy to the one without vortices. In that case, the existence of the pairs as the ground state will not depend 
only on the magnetization of the dipole but also on the parameters of the superconductor (i.e. A,^) which influence 
the equilibrium phase diagram via the self-energy of the vortices. 

From Fig. ^(b) it is clear that in the case of Ly — Lay = 1 there exists a critical MD magnetization for which the 
annihilation barrier disappears. More generally, the critical value of the magnetic moment of the MD, in the thin SC 
film limit, follows from Eq. ( p^ (by the condition that the function has no solution) and is given by 

m* = 1.08529 ^"^""^ . (14) 

Ly ~\~ Lay 

From this equation the critical values of the vorticity, when attraction becomes stronger than the pinning, can also 
be estimated. For Ly — Lav we have m* — 1.08529Lyl/2 and using the fact that <i>+/<i>o — 2m/l for an in-plane MD 
we obtain the critical condition $+/$o — l-08529Lt,, where denotes the flux through the region of the positive 
stray field of the dipole at the SC surface. This implies that there is a connection between the appearance of stable 
vortex-antivortex pairs and the amount of penetrating flux of thn magnetic dipole. Notice that the critical is 
not exactly quantized in units of <[>o which is a mesoscopic effect.E3 The quantization condition ^ = <f^ A ■ dl = L$o 
cannot be used because it is not possible to construct a contour C around the positive stray field region where the 
current is zero (see Fig. ^(b)). The inhomogeneous magnetic field of the dipole stimulates the creation of a vortex 
and an antivortex at opposite poles. These vortices cannot be treated independently and it turns out that larger flux 
is needed in order not only to create this pair but also to keep them apart. 

As shown above, for the case of the experimentally important thin SC film approximation, the position of the 
vortex-antivortex pair with respect to the dipole can be obtained from Eq. ([l^). This position depends both on 
the vertical position of the dipole and its magnetization strength. For fixed magnetic moment of the dipole, the 



vortex-antivortex pair can only be stabilized for I below some critical value which is obtained from Eq. (14). If we put 
TO = TO* in Eq. ([T^), we obtain the simple dependence p* = 0.4904/*. The dependence of p* on I is shown in Fig. ||, 
for different values of the magnetization, where also the stability region of the vortex-antivortex pair is indicated. One 
can see three regions in this diagram: i) the region where the vortex-antivortex pair is unstable, which is bounded by 
p* = 0.4904Z, for any value of the magnetization of the MD, ii) the region of stability, and iii) the forbidden area, for 
p* > 1.272L The latter condition follows from Sec. Ill (see Fig. ^(b)). Notice that with increasing distance between 
the dipole and the SC film, the pinning sites move further from the center of the dipole, up to a certain point when 
the interaction between vortices overwhelms the pinning force. After that, the distance between vortex and antivortex 
decreases and they finally annihilate. The maximum distance p* follows from Eq. ( p^ ) by taking the derivative with 
respect to I and leads to p* = 0.7071lmax (open dots in Fig. H), which again corresponds to the single vortex 
situation for I » d (see Fig. F 



B. Out-of-plane magnetized dipole 

Due to the strong field inhomogeneity and the reversal of the direction of the magnetic field in the vicinity of the 
MD with out-of plane magnetization, and the fact that the net magnetic flux due to the MD in the SC plane equals 
zerO|^ne expects that vortex-antivortex conflgurations might be stable in such applied field, as predicted earlier in 
Ref.u for a finite size superconductor. As shown in Fig. |^, an individual vortex is strongly attracted by the MD and 
the antivortices are repelled. Therefore we artificially put one vortex with vorticity Ly under the MD and assume 
the existence of a ring of Uay single antivortices around it, with radius Pt,.E3 Adding the interaction between each two 
vortices to our previous expression for the interaction energy, we obtain the total interaction energy 
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where the first two terms describe the dipole-vortex and dipole-antivortex interactions, the third term is the vortex- 
antivortex attraction term, and the last term is the repulsion energy between antivortices. The function Fy_av is 
given by Eq. (|2|). Although the vortex is attracted by the antivortices, annihilation is prevented by the repulsion 
between the antivortex and the MD (for example, see Fig. ^ for Ly = n^, = 1). Naturally, this energy barrier becomes 
smaller with decreasing magnetic moment of the dipole. Nevertheless, if a barrier exists, the interaction energy shows 
one local minimum, meaning that antivortices would not be repelled to infinity but to a certain point. The position 
of this local minimum we obtained in the same manner as previously, which within the thin film approximation, leads 
to the equation 

^""""^i, - 2Ly - nay + 1. (16) 

From this expression, the threshold value of the magnetic moment (when the potential barrier appears) is obtained 
as 



For lower values of the magnetic moment, the energetic barrier between vortices disappears and annihilation can not 
be prevented. Analogously, for fixed magnetic moment, increasing vorticity Ly of the vortex will make the attraction 
stronger and for a certain value of Ly, antivortices will be able to overcome the barrier. This critical value of vorticity 
can also be estimated from Eq. (p^. One should notice that we leave the possibility of Ly ^ Uay, which corresponds 
to the experimental situation when first positive external flux lines are pinned by the magnetic center, and then the 
polarity of the applied field is changed. Using <I>+/<i>o = 4\/37r/9 m/l for out-of-plane magnetic dipole polarization, 
we find for m — m* that the critical condition becomes <i>+/$o = (2i„ — Uay + 1) /2, where denotes the fiux 
through the region of positive stray field of the dipole. One should notice that the first stable vortex-antivortex pair 
[Ly = Uay — 1) appears for <!>+ — $o, and further increase of vorticity is a quantized process, with A$+ = <I)q/2 
necessary for stability for each one-unit-increase of Ly and Uay Also in this case we cannot define a path around 
the positive stray field region where the superconducting current is zero (see Fig. ^(a)) and consequently the flux 
quantization condition does not apply here. 

In this treatment, we assumed the presence of a giant vortex under the dipole as an energetically preferable state. 
It is well known that in infinite superconductors, thus in the absence of boundaries imposing the symmetry of the 
superconducting state, the giant vortex splits into multivortices. Since in our case an inhomogeneous applied magnetic 
field dictates the behavior of superconducting electrons, it is not clear which state carries less energy. Which central 
vortex configuration is realized depends on the parameters of the superconductor, i.e. ^, A and d, which come into 
play through the self energy of the involved vortices. Therefore, for a particular superconducting film, we extend our 
approach to the case of multivortices surrounded by an multi-antivortex ring. Namely, we investigate the stability 
of N vortex-antivortex pairs symmetrically arranged around the dipole, where vortices sit on a ring with radius py, 
and antivortices occupy the corresponding positions on the ring with radius pay We apply the same approach as 
before, calculating the interaction energy in a similar manner as in Eq. (p^), where the interaction between each two 
vortices is included. After minimization of the interaction energy with respect to the parameters py and pay, we 
obtain numerically the energetically favorable positions of the vortices. The results are shown in Fig. |l^, for the case 
of iV = 3 (see inset of Fig. |l^. Increasing magnetization of the dipole increases the distance between vortices and 
antivortices in a way that vortices come closer to the dipole and antivortices are repelled further awa^ It should be 
noted that a certain critical value of the magnetization is needed to prevent annihilation. In Fig. |l^ we give this 
threshold value as function of I for different values of N . One should note that this value of magnetic moment implies 
again the quantization of the penetrating fiux (as in the case of a giant vortex surrounded by antivortices), but with 
A$+ — 1.0489<I'o necessary for stability of the first pair, and slightly decreasing for additional pairs. For > 4, A$+ 
becomes smaller than the fiux quantum. The interpretation of this behavior can be that multivortices eventually join 
into a giant vortex, or that a different geometry of the superconducting state appears. For example, for larger number 
of pairs, due to the large linear density of vortices along the ring, they could rearrange, forming more than a single 
ring. 

Fig. ^ gives only the critical condition for stability of vortex-antivortex pairs. However, which configuration has 
the lowest energy and is thus the energetically preferable state cannot be inferred from this figure. In order to compare 




(17) 
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the energies of states with different number of pairs, i.e. N, we include the self-energy of individual vortices in the 
calculation 



$2 



167r2A 



d In— + 2tanh ^— 



In { 1 + cotanh [ — 



(18) 



The numerical results, as referred to the Meissner state, are shown in Fig. for a thin SC film with A/^ = 10. 
With increasing magnetic moment of the dipole the energetically favorable state goes through successive states in 
which N increases. If we calculate again the flux A$+ which now corresponds to the appearance of the next N 
state as the ground state, we find that A<i>+ = 1.97$oj for the appearance of the first pair, and for N > 2, the 
additional flux slightly decreases with increasing N, starting from A$+ = 1.089$o, for N = 2. Larger flux needed 
for the first vortex-antivortex state can be explained by the fact that this cylindrically asymmetrical state appears in 
a symmetrical magnetic potential. For N > 2, the results correspond to our previous analysis. It should be stressed 
that these results depend on the parameters of the superconductor, namely A and ^, which come into play through 
the self-energy of the vortices. This is in contrast to the interaction energy (see the inset of Fig. ^2|) which, in the 
London approximation, is independent of ^. 



V. CONCLUSION 



To summarize, we applied the London theory to investigate flux pinning in SC films due to the presence of a 
magnetic dipole situated above (or under) the SC. Depending on the direction of the dipole magnetic moment, we 
obtained exact analytic expressions for the MD-vortcx interaction energy and screening currents. We obtained the 
asymptotic behavior of the interaction potential and the induced currents for specific values of the involved parameters. 
We calculated the pinning potential for both an in- and out-of-plane magnetized dipole. Our results show that an 
out-of-plane magnetized dipole attracts a vortex if aligned parallel to it, and opposite, for anti-parallel alignment the 
vortex is repelled. This is a consequence of the mono-directional superconducting current induced in the SC for out- 
of-plane magnetization of the MD. However, for in-plane magnetization, the dipole-vortex interaction shows a dual 
behavior, namely, the vortex is attracted to the negative pole of the MD and repelled on the other side. Moreover, 
the position of the pinning site depends on the position of the MD and thickness of the SC. We calculated these 
dependences and showed that it is linear for thin superconductors, or large MD-vortex distances. Due to the dual 
behavior of the pinning potential, we explored the possible co-existence of vortices and antivortices in such systems. 
The total interaction energy calculation leads to the conclusion that the vortex and antivortex are separated by an 
energy barrier due to the short range interaction with the dipole, and therefore, these pairs could be stable. Both in- 
and out-of-plane magnetized dipoles are able to keep these vortices apart. We calculated analytically the interaction 
potential in the presence of vortex-antivortex pairs (or giant vortex-single antivortices) and gave estimates of the 
parameters necessary for stability of such fascinating configurations in a thin SC film. 
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APPENDIX: THE MAGNETIC FIELD OF A MAGNETIC DIPOLE IN THE PRESENCE OF A SC FILM 



The magnetic field of a vortex hy = rotA^, which is perpendicular to the plane of the film, is determined as the 
solution of the system of equations 

A^{p, z) + Tot rot Ay {p,z) ^ ^ {p- p^,) , \z\ < d/2, (Ala) 
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V2A„(p,z) = \z\>d/2, 



(Alb) 



with the following boundary conditions: i) the continuity of the vector potential components at z — ±(i/2, and ii) 
their vanishing far from the superconductor (at |^| — > oo). Here ^ (p — p^) takes into account the vortex, where py 



is the vortex position in the film plane 
rot<l = L$o<5(p - Pv)) 



In cylindrical coordinates "I>p 
The solutions of Eqs. (|AlaD-([Arb[) arJ^ 



0, $^ = L$o/27r |p - (so that 



42 (p, 



2ttX 



dq 



Jl (qR/X) 



kq cosh (kz/X) 



>sinh(M/2A) 



(A2a) 



2ttX 



Jl jqR/X) 
dq — cxp 



2\z\ 



2X 



(A2b) 



Here k = (l + g^)^ ^ , Q — k[k + gcoth {kd/2X)] , R = |p — py\ = [p"^ + p1 — '^PPv cos {ip — ip^)] ^^'^ , (a;) is the 
Bessel function, index "i(o\L. denotes the field inside (outside) the superconductor. The components of the vortex 
magnetic field are given by o 



1/2 



2ttX^ 



Ko 



dq q^ 
~kQ 



Jo 



qR 
T 



cosh {kz/X) 
sinh {kd/2X) 



(A3a) 



^ L<^>o f dqq f qR\ sinh(fcz/A) 
hil {P, - 2^ y -Q-J^ [- ) 3inh(fcrf/2A) ' 







L$0 [dqq (qR\ ( 2\z\~d \ 
h[} ip, ^)-^J -qJo [- ) exp [-q^j^ ) , (A3c) 







(P, ^) - ^ sgn(.) / ^ Jl ( ^) exp f-.^) , (A3d) 



2ttX^ " ' 'J Q ^\X ; ^ V 2A 



hi% = /il;^ = 0, (A3e) 

where Kq (x) is the MacDonald function. 

The magnetic field of a dipole, situated at the-jdistance / above the film, can be written in the following way (here 
and below I — \a\ ~ d/2) (see, for example, Ref.EII) 

{hd [p, z) + (p, z) , z> d/2, 
hi'} {p, z) , -d/2 <z< d/2, (A4) 

ht]-{P,z), z<-d/2, 

where hd (p, z) is the direct contribution from the dipole, h!"^ ^ (p, z) is the induced field above the superconductor, 
hm (p) is the dipole field penetrating inside the superconductor and (p) denotes the field of the dipole below the 
film. The vector potential Ad (p, z) is defined by the equation 

rot Ad (p, z) — iiTrnS (p) S {z — a) , (A5) 



where rh is the magnetic moment of the dipole. The potential Am (p) is found as the solution of Eq. ( Ala ) with zero 
right side and (p) as solution of Eq. (Alb). The integration constants have to be found through similar boundary 
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conditions as for the vortex field distribution: i) the continuity of the vector potential components at z = ±d/2 and 
ii) their vanishing at \z\ oo. 

Wc consider two orientations of the MD magnetization. 

I) The dipole is directed perpendicular to the film plane, i.e. nip = =0, rUz = m. In this case the vector 
potential has only an azimuthal component and is described by the following expressions: 



Adv{p,z) = ^ dq qJi (^y) exp (^-j \z - a\j = - 



mp 



p' + iz- ay 



3/2 ' 



(A6a) 



A%{p,z) = ^l^^dqqJ,{^) 



Bi (g) exp 



kz\ ^ , X { kz 
- +i?r(.)exp - 



(A6b) 



^m^,±(P> J2 '^1 (y) ^2,± (g) exp (-y) 



^m^ (P> ^) = ^3 / ^^"^^ (t) 



Bi (q) exp ( -y j + Di (q) exp 



(A6c) 



where Bi{q) = S(g)g (fc - g) exp (-/crf/2A) , Di{q) = B{q)q {k + q) exp {kd/2X) , B2.+ (q) 
sinh (fcrf/A) exp (gd/2A) , and ^2 _ (g) = 2B{q)kqexp {qd/2X) , with B (q) 

2m exp {—ql/X) / ^{k + qf exp (kd/X) — {k — q)^ exp {—kd/X)j . 

Consequently, we obtain the following components of the dipole magnetic field 



(A7a) 



h%{p,z) = ^ dqkqJ,[^) 



Bi {q) exp ( -y j - Di (q) exp f y 



(A7b) 



hdz{p,z) = ^ g^Jo (y) exp |z - a| 



2{z- ay - p2 



p2 + - ay 



5/2 ' 



(A7c) 



hdp {p, z) = 



msga.{z-a) f°° ^ (qp\ / q 



3mp {z — a) 



p' + iz- ay 



5/2 •■ 



(A7d) 



(P> ^) = ^3 / (t) ^^'^^ ^''P (~t) ' 



(A7e) 



C± (z^' = ^ X ^'-^i (f) ^2,± (a) exp (-^) , 



(A7f) 



/,(!) = /,(2) 



(A7g) 



II) The dipole is magnetized in the direction parallel to the SC film plane (in-plane magnetization). Following the 
same procedure as above, we obtain 



7 / ^ mcostp [°° ^ 2 / 9i A r T (<IP\ T ('iP\ 
hdp{p,V>,z) = dq q exp (^-- \z~a\j " "^o (,yj 



2p2 -(z- ay 



cos ip 



p^ + {z- ay 



5/2 



(A8a) 
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msmip 



dq gcxp {^-'^ \z - a|j Ji 



msmip 



p'^ + {z — ay 



3/2 ' 



(A8b) 



m cos Lp 
A3 



3mp |z — a| cos ip 



+ {z - ay 



5/2 ' 



(A8c) 



h% (p, ^, z) = ^ kq J2 ) - Jo (y ) 



Z?! ((7)exp ( y ] - -Si (q)exp ^-y 



(A8d) 



/^m^ (P> z) = 5^ kJi (y ) 



Di (g)exp ( y ] - -Si (?) exp ^-y 



(A8e) 



n-lnz (P> ^1 z) ^ ^ / dqq Ji[ — 



Di (q) exp ( y j + Bi (q) cxp ^-y 



(A8f) 



'imp.i (P,'/',^) = / rfg g^cxp (-y) J2(y)-Jo(y) B2,±{q), 



(A8g) 



Simp 



y^ qexp (^-y^ Ji (^y^ B2,± (g) , 



(A8h) 



(A8i) 
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FIG. 1. An oblique view of the systems under investigation with (a) out- and (b) in-plane directed magnetic dipole above a 
superconducting film interacting with a vortex. 

FIG. 2. The interaction energy of the vortex with an out-of-plane magnetic dipole (MD) in units of Fo = $2/(7rA) as a 
function of the position of the vortex (pt,): a) for different position of the MD; b) for different thickness of the superconductor; 
and c) the contour plot of the interaction potential. Inset in (a) illustrates the component of the total MD magnetic field in the 
presence of the SC, tangential to the SC surface {ho = mo/A^). Inset in (b) depicts the asymptotic behavior for specific values 
of the parameters. 

FIG. 3. The interaction energy of the vortex with an in-plane magnetic dipole (MD) in units of Fo = as a function 

of the position of the vortex (p„): a) for different vertical position of the MD; b) for different thickness of the superconductor; 
and c) the contour plot of the interaction potential. Open circles denote the energetically favorable position of the vortex. 

FIG. 4. The position of the vortex (p*) as function of the thickness of the superconductor (d) and distance between the 
in-plane magnetized dipole and the SC surface {I = a — d/2). p* is scaled by / or by A (inset of (b)). 

FIG. 5. The screening current density in the SC film with d = 0.5A for: a) out-of-plane and b) in-plane magnetized magnetic 
dipole with m = mo positioned at o = l.OA above the superconducting film. 

FIG. 6. The total interaction energy (in units of Fo = 3>o/(7rA)) between the in-plane magnetic dipole (MD) and the vortex 
and the antivortex placed at opposite sides of the dipole for different values of the magnetic moment m. 

FIG. 7. The interaction energy of a vortex-antivortex pair and an in-plane MD, as function of the position of the vortex 

(pav = —pv), along the direction of the magnetization of the dipole for: a) different values of the vorticity of the vortex and the 
antivortex, and b) for single-quantized vortices but different value of the magnetization of the dipole. The open circles denote 
the local minimum of the interaction energy. 

FIG. 8. The position of the vortex-antivortex pair (p*) as function of the distance between the in-plane magnetized dipole 
and the SC surface {I = a — d/2), for different values of the magnetization of the dipole. 

FIG. 9. The total interaction energy between an out-of-plane magnetic dipole (MD), the vortex (-I-) under it and the 
antivortex (-) at a distance p„. This illustrates vortex-antivortex pair stabilization by the short range interaction with the MD. 
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FIG. 10. The energetically favorable position of the vortex (+) and the antivortex (-) in the case oi N = 3 pairs placed 
symmetrically around the out-of-plane magnetic dipole (MD). The inset illustrates the configuration. 

FIG. 11. The stability phase diagram for A'^ vortex-antivortex pairs around the out-of-plane magnetized dipole, as function 
of the distance between the dipole and the SC surface (Z) and the magnetic moment m. Solid lines denote the critical value of 
the moment needed for the stability of such vortex-antivortex configurations. 

FIG. 12. The free energy difference between the states with N vortex-antivortex pairs around the dipole and the Meissner 
state, for a thin SC film with A/^ = 10 and the out-of-plane dipole positioned at l/\ = 1.0. Inset shows the total interaction 
energy in the corresponding A/^-states. 
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Fig. 1 Milosevic et al. 




Fig. 3(c) Milosevic etal. 
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Fig. 4(b) Milosevic etal. 
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